Multifractal wave functions of simple quantum maps 
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We study numerically multifractal properties of two models of one-dimensional quantum maps, 
a map with pseudointegrable dynamics and intermediate spectral statistics, and a map with an 
Anderson-like transition recently implemented with cold atoms. Using extensive numerical simula- 
tions, we compute the multifractal exponents of quantum wave functions and study their properties, 
with the help of two different numerical methods used for classical multifractal systems (box-counting 
method and wavelet method). We compare the results of the two methods over a wide range of val- 
ues. We show that the wave functions of the Anderson map display a multifractal behavior similar 
to eigenfunctions of the three-dimensional Anderson transition but of a weaker type. Wave functions 
of the intermediate map share some common properties with eigenfunctions at the Anderson tran- 
sition (two sets of multifractal exponents, with similar asymptotic behavior) , but other properties 
are markedly different (large linear regime for multifractal exponents even for strong multifractality, 
different distributions of moments of wave functions, absence of symmetry of the exponents). Our 
results thus indicate that the intermediate map presents original properties, different from certain 
characteristics of the Anderson transition derived from the nonlinear sigma model. We also discuss 
the importance of finite-size effects. 

PACS numbers: 05.45.Df, 05.45.Mt, 71.30.+h, 05.40.-a 
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I. INTRODUCTION 

Multifractal behavior has been observed in a wide va- 
riety of physical systems, from turbulence [l| to the stock 
market [2[ and cloud images 0]. It has been recog- 
nized recently that such a behavior can also be visible 
in quantum wave functions of certain systems. In partic- 
ular, wave functions in the Anderson model of electrons 
in a disordered potential are multifractal at the metal- 
insulator transition (see e.g. 0t9)- Similar behaviors 
were seen in quantum Hall transitions [7| and in Random 
Matrix models such as the Power law Random Banded 
Matrix model (PRBM) @,|3 and ultrametric random ma- 
trices [lO| . Such properties have also been seen in diffrac- 
tive systems [111 ] and pseudointegrable models, for which 
there arc constants of motion, but where the dynamics 
takes place in surfaces more complicated than the invari- 
ant tori of integrable systems [l^l- In all these models, 
this behavior of wave functions came with a specific type 
of spectral statistics, intermediate between the Wigner 
distribution typical of chaotic systems and the Poisson 
distribution characteristic of integrable systems [12, [3 ■ 
Recently, a new model of one-dimensional quantum "in- 
termediate map" which displays multifractal behavior 
was proposed |14{ , and a version with random phases was 
shown semi-rigorously to display intermediate statistics 
[15|. This model is especially simple to handle numeri- 
cally and analytically, and displays different regimes of 
multifractality depending on a parameter [16|. In par- 
allel, new experiments allowed recently to observe the 



Anderson transition with cold atoms in an optical po- 
tential [13, KLSl u sing a one-dimensional "Anderson map" 
proposed in [l9| . 

Although much progress has been made in the study of 
these peculiar types of systems, several important ques- 
tions related to multifractality remain unanswered. In 
particular, many results were derived or conjectured in 
the framework of the Anderson model, and their applica- 
bility to other families of systems is not known. Also, the 
precise link between the multifractal properties of wave 
functions and the spectral statistics is not elucidated. 

In order to shed some light on these questions, we 
systematically investigate several properties of the wave 
functions of the intermediate quantum map of [14l llSj , 
and compare them to results obtained for the Anderson 
map. In these one-dimensional systems, wave functions 
of very large vector sizes can be obtained and averaged 
over many realizations. This enables to control the errors 
and evaluate the reliability of standard methods used in 
multifractal analysis. This also allows to check and dis- 
cuss several important conjectures put forth in the con- 
text of Anderson transitions. Our results also enable to 
study the road to asymptotic behavior in such models, 
giving hints on which quantities are more prone to finite- 
size effects or can be visible only with a very large number 
of random realizations. 

The manuscript is organized as follows. In Section II, 
we review the known facts and conjectures about quan- 
tum multifractal systems, which were mainly put forth in 
the context of the Anderson transition. In Section III, we 



present the models that wih be studied throughout the 
paper. Section IV discusses the numerical methods used 
in order to extract multifractal properties of the models. 
Section V presents the results of numerical simulations, 
allowing first to compare the different numerical methods 
of Section IV, and then to test the conjectures and re- 
sults developed in the context of the Anderson transition 
to the two families of models at hand. 



II. EARLIER RESULTS AND CONJECTURES 

We first recall some basic facts about multifractal anal- 
ysis. Localization properties of wave functions of compo- 
nents ipi, i = 1, ..., N, in a Hilbert space of dimension N, 
can be analyzed by means of their moments 
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The asymptotic behavior of the moments ([ij for large N 
is governed by a set of multifractal exponents Tq defined 
by Pq ex N^'^'' , or by the associated set of multifractal 
dimensions Dq = Tq/(q~l). Equivalently, the singularity 
spectrum /(a) characterizes the fractal dimensions of the 
set of points i where the weights jV'iP scale as N~". It 
is related to the multifractal exponents Tq by a Legendre 
transform 
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Compared to classical multifractal analysis, the quan- 
tum wave functions in Hilbert space of increasing dimen- 
sions are considered as the same distribution at smaller 
and smaller scales. This allows to define properly the 
multifractality of quantum wave functions, although at a 
given dimension they correspond to a finite vector. 

In many physical instances, only a single realization 
of the wave function is considered. However, as soon 
as several realizations are considered, as is the case in 
presence of disorder, the moments ([T|) are distributed ac- 
cording to a certain probability distribution, and mul- 
tifractal exponents depend on the way ensemble aver- 
ages are performed, and in particular on the treatment 
of the tail of the moments distribution. If the tail de- 
creases exponentially or even algebraically with a large 
exponent, the two averages should give the same answer. 
On the other hand, if the moments decrease according 
to a power- law with a small exponent, the average (Pq) 
will be dominated by rare wave functions with much 
larger moments (whose magnitude directly depends on 
the number of eigenvectors considered), while the quan- 
tity Pq^^ = exp(lnPq) will correspond to the typical 
value of the moment for the bulk of the wave functions 
considered. To each of these possible averaging proce- 
dures corresponds a set of multifractal exponents, defined 
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As soon as averages over several realizations arc made 
there can be a discrepancy between Tq and t'^'p. Histor- 
ically this effect was seen in the context of the Ander- 
son transition [6|, |2l| and was ver y re cently confirmed by 
the analytical calculations of [23, l2j| in the same model. 
More specifically, it was shown that the distribution of 
the normalized moments yq = Pq/Pq^^ is asymptotically 
independent of N and has a power-law tail 
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The multifractal exponents Tq 
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In the case of heavy tails Xq < 1, the averages (Pq) and 
exp(lnPq) yield different exponents. 

This phenomenon has a counterpart in the singular- 
ity spectra /(a) and /'yP(a). While /*-^P(a) cannot take 
values below and terminates at points a± such that 
/'yp(a_l-) = 0, the singularity spectrum /(a) continues 
below zero. The two spectra coincide over the inter- 
val [a+,a_]. It can be shown that outside the interval 
[(7_,g+] the set of exponents r'^P are given by the lin- 
ear relation t'-^p = qa^ for q > q+ and t'^p = qa^ for 
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, it was stated that the exponents Tq and t'^p can 
be related through the following relation which depends 
on the tail exponent of the moments distribution Xq: 
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Equation ([6]) was analytically proven for PRBM [8| for 
integer values of Xq and also in the limit of large bands 
for q > 1/2. It remains unclear to which extent Eq. ([6]) 
is valid for other types of systems. A consequence of the 
identity ^ ioi q > (7+ is that Xq = q+/q. In the regime 
of weak multifractality where Dq is a linear function, the 
identity ([6]) implies that Xq — [q^jqY for q^ < q < q+ 
[all] (see Subsection V.E. for details). 

Finally, a further relation that we wish to investigate 
in the present paper has been predicted based on the 
nonlinear sigma model and observed in the 3D Anderson 
model at criticality and several other systems. It is a 
symmetry relation of multifractal exponents [25| , which 
can be expressed as 
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with Aq = Tg — g + 1. For the singularity spectrum it 
gives /(2 -a)= f{a) + l-a. 

The validity of these different relations will be inves- 
tigated on two particularly simple models of quantum 
maps that we describe in the next section. 



III. MODELS 



A. Intermediate map 



The system is a generalization of the quantum kicked ro- 
tator model, and is described by a unitary operator which 
evolves the system over one time interval: 



The properties of section |TT] have been first observed 
for wave functions in the 3D Anderson model at critical- 
ity. In the present paper the first model we will consider 
is a quantum map whose eigenfunctions display similar 
multifractal properties in momentum representation. It 
corresponds to a quantization of a classical map, defined 
on the torus by 
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p ~ p + ^ (mod 1) ; q ^ q + 2p (mod 1) 



(8) 



where p is the momentum variable and q the angle vari- 
able, while p and q are the same quantities after one iter- 
ation of the map. The corresponding quantum evolution 
can be expressed as 
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in operator notation, or equivalently as a, N x N matrix 
in momentum representation: 
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(10) 
with < fc, fc' < A^ — 1 [ll|. For generic irrational 7, the 
spectral statistics of U are expected to follow Random 
Matrix Theory. When 7 is a rational number, 7 = a/b 
with a, b integers, spectral statistics are of intermediate 
type and eigenvectors of the map display multifractal 
properties. In order to study the effect of ensemble aver- 
aging on multifractal exponents, we will consider a ran- 
dom version introduced in [15[, where the phases 27rfc^/7V 
are replaced by independent random phases 0^. We will 
also present numerical results obtained for the initial map 
given by (fTOj) . 



B. Anderson map 

An important system where multifractal wave func- 
tions have been observed corresponds to electrons in a 
disordered potential in dimension three. Indeed, the An- 
derson model which describes such a situation displays 
a transition between a localized phase (exponentially lo- 
calized wave functions) and a diffusive phase (ergodic 
wave functions) for a critical strength of disorder. At 
the transition point, the wave functions display multi- 
fractal properties Q and the spectral statistics are of 
the intermediate type [26|. In order to compare this 
type of system with the previous one, we have studied 
a one-dimensional system with incommensurate frequen- 
cies, which has been shown to display an Anderson-like 
transition [l^. In [T^ it was shown that it can be im- 
plemented with cold atoms in an optical lattice, which 
enables to probe experimentally the Anderson transition. 



with V{9,t) — fc(l 4- ecoscJitcosw2i) cos^ (here time t 
corresponds to number of kicks) . Here lo\ and 102 should 
be two frequencies mutually incommensurate. Following 
|17| we chose in the simulations e = 0.75, wi = 27rA~^ and 
L02 = 27rA~^, where A = 1.3247 ... is the real root of the 
cubic equation cc"^ — .t — 1 = 0. In |la| it was shown that 
this system displays an Anderson-like transition at the 
critical value kc ~ 1.81, but multifractality of this system 
was up to now not verified. The function Hq^h) can be 
chosen either by taking Ho{n) = n^/2 (free evolution) 
or as in the preceding case by replacing this quantity 
by independent random phases uniformly distributed in 
the interval [0, 2tt[. This is actually what we chose to 
do in our numerical simulations, in order to increase the 
stability and accuracy of the numerical results. 

As a wave packet spreads slowly at the transition, one 
has to iterate the map for a long time in order to obtain 
data on a sufficiently large wave function. We found that 
in order to reach vector sizes of order 2^^, it was necessary 
to iterate the map up to t = 10^. Such values are not 
realistic for experiments with cold atoms (limited to a few 
hundreds kicks) but allow to obtain more precise results. 



C. Variations on the models 

In both models the evolution of the system during one 
time step has the form of an operator diagonal in momen- 
tum (kinetic term) times an operator diagonal in position 
(kick term). In both cases, it has been common practice 
in the field to replace the kinetic term by random phases. 
This allows to obtain a similar dynamics but with a more 
generic behavior. Moreover, it enables averages over ran- 
dom phases to be performed, which makes numerical and 
analytical studies more precise. In contrast, in experi- 
ments with cold atoms it is easier to perform iterations 
with a true kinetic term rather than with random phases. 
In order to assess the effect of this modification, we will 
therefore use both approaches in the study of multifractal 
properties of wave functions. 

Many works on multifractal wave functions have fo- 
cused on eigenstates of the Hamiltonian (see e.g. Q)- For 
the intermediate map, the evolution operator has eigen- 
vectors which can be numerically found and explored. It 
is also possible to evolve wave packets, for example ini- 
tially concentrated on one momentum state, and to study 
the multifractality of the wave packet as time increases. 
This corresponds more closely to what can be explored 
in experiments. For the intermediate map, this process 
can be understood as the dynamics of a superposition of 
eigenvectors of the evolution operator. However, in the 
case of the Anderson map, the evolution operator is time 



dependent (the continuous time problem being not pe- 
riodic), and the connection with eigenvectors is lacking. 
In the following, we will explore and compare the multi- 
fractality of both eigenfunctions and time-evolved wave 
packets for the intermediate map. 



IV. NUMERICAL METHODS 

As is well-known, the numerical estimation of multi- 
fractal dimensions is very sensitive to finite-size effects. 
In the present work we have analyzed and compared dif- 
ferent numerical methods in order to compute accurately 
the multifractal exponents. In this section we briefly re- 
view the methods we used. 



A. Box-counting method 

The most straightforward method is to compute di- 
rectly the moments of the wave functions through the 
scaling of the moments ([T]) given by Pq = X^i iV'iP'^- If 
the scaling ^ holds true, then log(Pq) is a linear func- 
tion of log N and its slopes yield the exponents Tq . For 
g < 0, coarse-graining over neighboring sites is necessary 
in order to avoid instabilities due to very small values of 

A variation of the moment method is the box-counting 
method. It consists in taking a vector of fixed size N 
and summing the ji/'iP over boxes of increasing length. 
If iV = 2", then we define 



q,k 




(12) 



which corresponds to averaging the measure over boxes 
of length 2'^ , < k < n. Starting from fc = 2 allows to 
smooth out exceedingly tiny values of IV'iP which other- 
wise yield non-accurate estimates of Dq for negative q. 
The two methods above give similar results, as was also 
carefully checked for the Anderson model in |20| . We will 
therefore in the following present results using the box- 
counting method to assess the properties of this type of 
procedures. 



B. Wavelet transform method 

Recently, alternative procedures to compute the mul- 
tifractal spectrum based on the wavelet transform were 
developed [23| . The wavelets form a basis of functions as 
does the Fourier basis and a function expressed in this 
new basis gives the wavelet transform (WT). Unlike the 
Fourier basis wavelets are localized both in position and 
in momentum space (or time and frequency space) . They 
are therefore suitable to probe the local variations of a 
function at different scales. They have become an essen- 
tial tool for image and sound processing and compression. 



A wavelet basis is constructed from a single function g 
called analyzing or mother wavelet. The rest of the basis 
is constructed by translations and compressions (expan- 
sions) of the analyzing wavelet g. The translations define 
a space variable while the compressions define the scale 
at which the function is analyzed. We define the WT of 
a (real) function h as 



Th{A,B) 



dx h{x) g 
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(13) 



where A is the scale variable and B is the space variable. 
As a consequence, Th{A^ B) can be interpreted as a mea- 
sure of how close the function h is to the mother wavelet 
at point B and at scale A. 

The Tq can be extracted from the WT in the following 
way. We define the partition function 



Z(g,A) = ^|r,.(A,BO 



(14) 



where (A, Bi)i are the local maxima at scale size A and q 
is real. It can be shown that Tq appears as the exponent 
in the power law behavior of 2(g, A) 
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This is essentially the method known as wavelet trans- 
form modulus maxima (WTMM) [13 • This method is 
developed for continuous wavelet function. 

If the function h is sampled as an A^-dimensional vector 
with N = 2" the wavelet transform can be discretized 
and implemented efficiently by a hierarchical algorithm 
[28| resulting in a fast wavelet transform, (FWT). The 
scale and space parameters take the values 
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(16) 
(17) 



respectively. Starting from a wave function ip and us- 
ing a proper normalization at each scale, we redefine the 
partition function as 



Z{q,A)=Y^ 



\T\.^\2{A,B,)\ 



(18) 



where again {A, Bi)i are the local maxima at scale A. As 
in the continuous case Z{q, A) exhibits the same power 
law behavior as Eq. (|15p . In the following, we will present 
results using this implementation of the wavelet method, 
using the Daubechies 4 mother wavelet [29| . 

We note that, although the partition function ([T8|) is 
the most standard, recently it was numerically observed 
[30| that for a complex multifractal wave function -0, the 
exponents Tq can also be obtained from the power law 
behavior of a modified partition function defined from 
the complex FWT of V' a-nd using 2q as exponents in 

Eq. m 



V. RESULTS 

A. Numerical computation of the multifractal 
exponents 
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FIG. 1: (Color online) Top panel: Instance of an eigenvector 
of the intermediate map with random phases for A'^ = 2^* 
and 7 = 1/3. Bottom panel: Instance of an iterate of the 
Anderson map with random phases for k — 1.81, A'^ — 2^^ 
and t = 10*; j?/'(0)) = \i) with i = 1024. 



In this subsection, we present numerical results corre- 
sponding to the multifractal exponents for the interme- 
diate quantum map model (fTO)) and the Anderson map 
([Tl]) . Examples of wave functions for both models are 
shown in Fig. [1] The two sets of multifractal dimen- 
sions Dq and D^J'^ were computed from log2(Pg,fc) and 
(log2 Pq,k) respectively for the box-counting method, and 
from log2 {Z[q, A)) and (log2 Z{q, A)) respectively for the 
wavelet method. For the Anderson map, as less realiza- 
tions of the random phases could be computed, the same 
set of random realizations were used for the two methods 
investigated in order to ensure that comparable quan- 
tities were plotted. Figure [5] illustrates the scaling of 
these quantities for the two models considered. It dis- 
plays log2(P5,fc) ^^ ^ function of the logarithm of the box 
size (top) and \og2{Z{q,A)) as a function of the loga- 
rithm of the scale parameter log2 A (bottom) for different 
values of q. We chose to show the scaling of these mo- 
ments since it is the worst case, the curves for (log2 Pq.k) 
and (log2 Z{q, A)) being always closer to linear functions. 
Nevertheless, Fig. [5] shows that the scaling is indeed lin- 
ear over several orders of magnitude. The slopes of the 
linear fits give the multifractal exponents. In the two 
methods, there is a certain freedom in determining the 
range of box sizes (or scales for the wavelet method) over 
which the linear fit is made (the one we chose is indicated 
by the shaded area in Fig. ^ . Usually for moderate val- 
ues of q in absolute value, the result does not depend 
very much on this choice. As can be seen in Fig. [21 the 
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FIG. 2: (Color online) Top panels: Average moments 
log2{P5,fc) as a function of k, the logarithm of the box sizes. 
Bottom panels: Partition function \og2{Z{q, A)) as a function 
of the logarithm of the scale A. The left panels correspond 
to eigenvectors of the Intermediate map with random phases, 
N = 2^*, and 7 — 1/3. Here the average is taken over 98304 
vectors (resp. 32768 vectors for the partition function). The 
right panels correspond to iterates of the Anderson map for 
k = 1.81, N = 2^\ IV'(O)) = \i) with i = 1024 and t = 10^ 
where the average is taken over 1302 vectors. The values cho- 
sen for q are : q = —2 (red diamonds), q — 2 (blue circles), 
and g = 6 (green squares). The gray shaded regions show the 
fitting interval. 



data are well fitted by a linear function in the range cho- 
sen; however, there is still an uncertainty on the slope, 
which usually gets larger for large negative q. In the next 
figures of this subsection, the uncertainty of the linear fit 
for the set of points chosen is plotted together with the 
mean value, in order to give an estimate of the reliability 
of the values obtained. 

The values of Dq and D^^p as a function of q are pre- 
sented in Fig. [3] for the random intermediate map and in 
Fig. [4] for the Anderson map. In both cases, we observe 
a spectrum typical of multifractal wave functions. In the 
intermediate case, we observe that the two methods give 
comparable results with small uncertainty, although it 
gets larger for large negative q. For the Anderson map, 
the uncertainty gets very large for q < —2, and besides 
it begins to depend strongly on the range of box sizes 
(box-counting method) or scales (wavelet method) cho- 
sen (data not shown). We attribute the larger uncer- 
tainty for the wavelet method to the seemingly stronger 
sensitivity of this method to exceptionally small values. 
We note that for the Anderson map only a vector size of 
2^^ could be reached numerically because of the compu- 
tational power required to iterate the Anderson map for 
long times. The discrepancy between the two methods 
is smaller for -D'^p, which can be understood by the fact 
that {Pq,k) or {Z{q,A)) are more sensitive to rare events 
than (log2Pq,fe) and (log2 2^(g, A)). 
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FIG. 3: (Color online) Top panel: Multifractal exponents Dq 
for eigenvectors of the Intermediate map with random phases, 
A'^ = 2^*, and 7 = 1/3. Empty/filled red circles correspond to 
the method of moments/wavelets. The (light blue, light red) 
shaded regions indicate standard error in the least squares 
fitting. The multifractal analysis was done over 98304 vectors 
with box sizes ranging from 16 to 1024 and scales ranging 
from 2 — 12 to 2 — 5. Bottom panel: Typical multifractal 
exponents D^J^ for the same data. In both panels, the gray 
solid line is the linear approximation 1 — g/3. 



Our data nevertheless show that the iterates of the 
Anderson map display multifractal behavior. In view of 
the recent implementation of such maps with cold atoms 
[17| , this indicates that in principle one can observe mul- 
tifractality of this map in cold atom experiments; we note 
that in [l3l properties of the wave functions of this ex- 
perimental system were investigated, but focused on the 
envelope of the wave packet. 

In both the intermediate map and the Anderson map, 
D*yp goes to a constant a+ for large q, which corresponds 
to the fact already mentioned that t^^ is expected to 
behave as r*^P = qa+ for q > q+. This will be studied in 
more detail in the next section. 

In jl6| we observed the existence of a linear regime 
around q = 0, with slope —1/6 for the random interme- 
diate map with parameter 7 = 1/6. This regime seemed 
to be valid in quite a large interval around zero. The 
data displayed in Fig. [3] obtained by two different large- 
scale computations confirm this result, and show that 
this regime exists for both types of averages. We will 



FIG. 4: (Color online) Multifractal exponents Dq and D'^^ 
for iterates of the Anderson map with k = 1.81, A'^ = 2 ^. 
Same convention as in Fig. [S] The multifractal analysis was 
done over 1302 realizations of size A^ — 2^^ with box sizes 
ranging from 4 to 512 and scales ranging from 2 — 9 to 2 — 4. 



come back to this property in section FVEI 

The data presented in Fig. |4] show that the Anderson 
vectors are less fractal than the eigenvectors of the in- 
termediate map for 6 = 3. This can be expected from 
the fact that Anderson vectors correspond to iterates of 
a quantum rnap, which in general are less fractal than 
eigenvectors |3l| . In order to compare similar quantities, 
we display in Fig. [5] multifractal dimensions for vectors 
obtained by iteration of the intermediate map for 6 = 3. 
One sees that using iterates instead of eigenvectors also 
in this case reduces the overall multifractality at a given 
q. Although a simple relation between multifractal ex- 
ponents of eigenvectors and iterates is still lacking, this 
indicates that it might be possible to use experimental 
results from cold atom experiments to infer properties of 
the eigenvectors of the 3D Anderson transition. 

In order to assess the effect of random phases, we 
present in Fig. [6] the results of numerical computation 
of Dq and D*^p for the intermediate map (fTO|) without 
random phases. Although spectral statistics for random 
and non-random vectors are very close [l^ , the obtained 
fractal dimensions are quite different. Moreover it can 
be seen that there still exists a difference between Dq 
and -D'-^P; although the map is not random any more. 
The discrepancy in the two sets of exponents was to our 
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FIG. 5: (Color online) Top: Multifractal exponents Dq and 
Z)*yp for the Anderson map with k = 1.81, A^ = 2", 
1 1/1(0)) = \i) with i = 1024 and t = 10*. Same convention 
as in Fig. [3l The multifractal analysis was done over 1302 
realizations with box sizes ranging from 4 to 512. Bottom: 
Multifractal exponents Dq and D^^ for iterates of the inter- 
mediate map with b = 3, iV = 2^^ |?/>(0)) = |i) with i = 2048 
and t = 10*. The multifractal analysis was done over 2000 
realizations with box sizes ranging from 4 to 512. In both 
cases we used the box-counting method. 



knowledge up to date observed only in disordered sys- 
tems, where rare events are created by specific realiza- 
tions of disorder. It is interesting to see that this effect 
can be observed in a dynamical system without any dis- 
order whatsoever. This discrepancy is due to the fact 
that the average in Eqs. (I3|)-(l4|) is performed over sev- 
eral eigenvectors of a single realization of the map, which 
gives a certain dispersion of the moment distribution. 
The average over eigenvectors in the intermediate map 
thus sufHces to create the effect even in a determinis- 
tic map. Therefore one can also observe the separation 
between the two sets of exponents in quantum systems 
without disorder. The rare events in this case correspond 
to rare eigenvectors of the evolution operator having large 
moments. 

The numerical results displayed in this Subsection in- 
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FIG. 6: (Color online) Top panel: Multifractal exponents Dq 
for eigenvectors of the intermediate map with (circles) and 
without (triangles) random phases. Bottom panel: Typical 
multifractal exponents D^'^ for eigenvectors of the interme- 
diate map with (squares) and without (diamonds) random 
phases. In both cases we used the box-counting method. Con- 
ventions and parameter values are the same as in Fig. O ex- 
cept that for non-random phases, a single realization of size 
N = 2^"^ was considered. In both panels, the gray solid line 
is the linear approximation 1 — q/3. 



dicate that mult ifr act ality is indeed present in all models 
considered. Furthermore, our results show that for an 
appropriate choice of a range of box sizes (box-counting 
method) or scales (wavelet method), both methods are 
based on curves well-fitted by linear function over a wide 
interval, with a small uncertainty on the exponent ex- 
tracted. We believe that the good agreement for q> —2 
between the two methods, and the small uncertainty 
found for the linear fit, indicates that our numerically 
extracted multifractal exponents are reliable (up to fi- 
nite size effects). For q < —2, the numerical uncertainty 
increases for decreasing q, and the two methods give in- 
creasingly different results. The results presented here 
show that our data arc still reliable, although less pre- 
cise, for eigenvectors of the intermediate map, even for 
negative q. However, in the case of the Anderson map, 
the uncertainty is too large to give reliable results for 
q < —2. This is certainly due to the fact that both the 
number of realizations and the vector sizes are smaller 
in this case, which makes it difficult to find reliable re- 
sults in the more demanding regime of large negative q. 



In the regime of g > —2, our results indicate that the 
two methods can be used and give similar results. In 
the more demanding cases (large negative q) we found 
the box-counting method more reliable and accurate, and 
therefore the numerical results presented in the following 
Subsections correspond to this method. 



B. Linear regime 

In [Ig] the first investigations of the multifractal ex- 
ponents for the random intermediate map showed the 
presence of a linear regime around q = over a relatively 
large range of q values. We recall that the map (fTO|) 
displays multifractal properties of rational values of the 
parameter 7 = a/b. Based on semi-heuristic arguments, 
this linear regime was shown to be described by: 
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The relation P^ enables to link the spectral statistics 
and the distribution of Dq around g = in a systematic 
way since both are controlled by the parameter b explic- 
itly. While one expects that some form of linear regime 
should exist over small intervals for any smooth curve, 
the extent of it in this particular model indicated a small 
second derivative near q = 0. This feature was seen in 
the PRBM model p] but in the regime of weak multifrac- 
tality where the derivative a.t q ~ oi Dq is very close to 
zero. 

Figure [7] displays the extent of the linear regime for 
the intermediate map for three values of the parameter 
7 = 1/b. The data presented show that the linear regime 
is present in all three cases, although its extent seems to 
be larger for large b (weak multifractality) . This indicates 
that the linear regime is a robust feature of the random 
intermediate map. To explore more precisely the depen- 
dence of this regime on the value of b, the inset of Fig. [7] 
shows the separation point between the actual exponent 
and the linear value, defined by a constant relative dif- 
ference (set to 1%), for all values of b between 6 = 2 and 
b = 13. The data presented show that indeed the extent 
of the linear regime grows with b, although the precise 
law of this growth is difficult to specify. 

These results correspond to the random intermediate 
map, where the kinetic term is replaced by random phase. 
It is important to explore also the behavior of the deter- 
ministic intermediate map, where the kinetic term is kept 
as a function of momentum. In this case. Fig. [7| shows 
the presence of a much smaller linear regime. This indi- 
cates a strong difference between the random model and 
the deterministic one, and that a large linear regime is a 
property restricted to a certain kind of models. 

The data for the Anderson map with random phases, 
shown in Fig. [8l also show a linear regime, comparable 
with the random intermediate map. Again, as the data 
correspond to iterates of wave packets, not eigenvectors, 
the multifractality is weaker than for other simulations 
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FIG. 7: (Color online) Difference Dq — Dg'^ for eigenvectors of 
the intermediate map with random phases for 7 = 1/3 (blue 
solid curve), 7 = 1/7 (green dotted curve) and 7 — 1/11 
(red dash-dotted curve). Blue dashed curve shows the same 
difference for the intermediate map without random phases 
for 7 = 1/3. Inset: separation points between Dg and D^" 
(green dots), determined by {Dg - D|,'")/(D, + D'^'') = 0.01, 
for eigenvectors of the intermediate map with random phases 
for different values of 7 = 1/b. Other parameters as in Fig.[51 
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FIG. 8: (Color online) Difference Dg 
Anderson map, with D^" defined by Dg" 
parameters as in Fig. [l] 



Dg" for iterates of the 



l + gD;(0). Same 



of the Anderson model using eigenvectors [6| . This might 
explain why the linear regime that is visible seems larger 
than for Anderson transition eigenstates. 



C. Average vs typical multifractal exponents 

In disordered systems, the statistical distribution of 
the moments of the wave function is responsible for a 
discrepancy between the multifractal exponents Dg and 
-D'^P calculated respectively by averaging over the mo- 
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FIG. 9: (Color online) Top panel: Exponents r, (red solid 
lower curve) and r^^^ (blue solid upper curve) as a function 
of q for eigenvectors of the intermediate map with random 
phases, A^ = 2^*, and 7 = 1/3. Bottom panel: same but for 
7 = 1/11. In both panels, the black dashed lines are linear 
fits of Tg^P at large \q\. The arrows indicate the critical q 
value determined by a;^ = 1 (see Subsection V.D), equal to 
2.51 for 7 = 1/3 and to 4.30 for 7 = 1/11. The slopes of the 
black dashed lines are respectively given by a_ ~ 2.28 and 
Qf+ ~ 0.32 for 7 = 1/3 and by a- ~ 1.57 and q+ ~ 0.60 for 
7 = 1/11. 



mcnts themselves or over their logarithm. The two sets 
of exponents are expected to match only in some region 
q G [q-,q+]. Outside this range, r^^P should follow a 
linear behavior. Figure [9] displays results for r, and t*^p 
for the intermediate map with parameter 7 = 1/& for two 
representative values of b. As expected from the PRBM 
model Q, in the case of weak multifractality (7 = 1/11, 
bottom panel of Fig. [9]), the range over which the expo- 
nents are equal is wider than for strong multifractality 
(7 ~ 1/3, top panel of Fig. [9|) . Beyond that interval the 
behavior of t'^p is linear for both values of 7. For positive 
q the linear tail appears around q ~ 2.5 for 7 = 1/3 and 
q ~ 4.3 for 7 = 1/11. The slopes of the linear tails give 
q;_ and a_|_. According to the theory (see Section II), 
these values of a correspond to the terminating point of 
the singularity spectrum /'yP(a) defined in Eq. ([2]). We 
obtained comparable results for the Anderson map (data 
not shown). 

In order to take a closer look at the discrepancy be- 
tween the two sets of exponents, we plot the difference 



FIG. 



10: (Color online) Top panel: Difference D, - 1)*^^ be- 



tween average and typical exponents for eigenvectors of the 
intermediate map with random phases for 7 = 1/3 (red solid 
curve) and 7 — 1/11 (blue dashed curve). Bottom panel: 
same figure for iterates of the Anderson map. Same parame- 
ters as in Fig. |4l 



between the exponents Dq and D^^p for both systems. 
It is clearly seen in Fig. [10] that in all cases the regime 
where the exponents coincide is only about —l<q<l. 
At this scale the separation between r^ and t'^p occurs 
around q = 1. In order to obtain more systematically the 
separation point, we have plotted in Fig. [Tl] the value of 
q defined by a constant relative value of the difference 
of the two exponents (set equal to 0.01) for intermediate 
maps with different parameters 7 = 1/6; this allows to 
get comparable data independently of the value of the 
exponents. The results show a clear linear scaling of the 
separation point with respect to b, obeying the formulas 
q-0.01 + 1 « -0.1(6 + 1) and q+0.01 - 1 ~ 0.15(6 - 1). 
Changing the threshold of relative value from 0.01 to 0.02 
gives also a linear scaling, with different slope (data not 
shown). 

The singularity spectrum /(a) defined in Eq. ^ is an 
alternative way to analyze multifractality and the dis- 
crepancy between the two sets of multifractal exponents. 
In Fig. [12] we show the singularity spectrum obtained 
for the intermediate map pop with random phases (in 
[iq similar curves were obtained using directly the box- 
counting method to compute /(a); here we use the Leg- 
endre transform of the exponents, obtaining similar re- 
sults). As expected /'^P(a;) terminates at points a^ and 
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FIG. 11: (Color online) Separation points g±o.oi between Dq 
and D*/P (green dots), determined by {Dg - Df^)/{Dq + 
D^q^) = ±0.01, for eigenvectors of the intermediate map with 
random phases for different values of 7 = 1/b. The blue solid 
lines are the linear fits g-0.01 ~ —0.106 — 1.10 and g+0.01 ~ 
0.156 + 0.86. Same parameters as in Fig. |4l 




FIG. 12: (Color online) Left panel: Singularity spectra /(a) 
and /*■*'''(«) for eigenvectors of the intermediate map with 
random phases for 7 = 1/3. The inset shows the linear be- 
havior of \/b/(a_ — Qf+) ~ 0.3076 — 0.037 as a function of 
6 (in 7 = 1/6). The values of a± have been extracted from 
linear fits of t^"^ at large |g|. Right panel: Singularity spec- 
tra /(a) and /*■*'''(«) for iterates of the Anderson map. In 
both panels, the singularity spectra are given in the range 
a{q = —16) < a < a{q — 16). 



a- given by the large-g slopes of t'^'p, while f{a) takes 
values below zero coming from statistically rare events. 
In [l6| , it was shown that a linear approximation for Dq 
yields a parabolic approximation for /(a), giving in turn 
a behavior a_— «+ ~ l/Vb. We checked this behavior for 
all values of b between b — 2 and b = 13, thus confirming 
the validity of this law, even beyond the linear regime 



(see inset of Fig. [T2|) . Figure [12] allows a more direct 
comparison between multifractality in the intermediate 
map and the Anderson map: the narrower /(a) curve 
for Anderson corresponds to a weaker multifractality. 



D. Moment distribution 

The discrepancy between the two sets of multifractal 
exponents observed in the previous section is due to the 
fact that the moments Pg defined by ^ have a statis- 
tical distribution with a certain width. In particular for 
multifractal measures the distribution of the normalized 
moments yq = Pq/Pq^^ is expected to have a power-law 

tail at large q as V{yq) ^ 1/j/g '; with an exponent Xq 
depending on q. In Fig. [T3] we show an example of the 
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FIG. 13: (Color online) Probability distribution of the loga- 
rithm of the moment Pq = P{q, 0) for q — 1.25 of eigenvectors 
of the random intermediate map for 7 — 1/3. Same parame- 
ters as in Fig. [3] The solid line shows a linear fit (in log scale) 
whose slope yields the tail exponent —Xq. 

distribution of the logarithm of the moments Pq for the 
random intermediate map. The distribution is indeed al- 
gebraic with a power law tail depending on q. While the 
linear behavior (in logarithmic scale) is clearly observed 
for small values of q (restricted to the range q > 1), this is 
not the case for larger q. We calculated the exponent Xq 
of the tail for a range of values of q where this exponent 
could be extracted. Results are displayed in Fig. [M] for 
the intermediate map and Fig. [15] for the Anderson map. 
The value of q where Xq ~ 1 should correspond to the 
value q+ where the two curves of Dq and ZJ'^p (or Tq and 
T*-^P) separate. As one can observe in Fig.[T4]that value of 
q is rather difficult to estimate numerically with sufficient 
accuracy as the exponents Xq did not converge to a defi- 
nite value at the largest vector size available {N = 2^'^). 
However the curves seem to yield an exponent equal to 1 
around q « 2.51 for 7 = 1/3 and 4.3 for 7 = 1/11. These 
values are indicated with black arrows in Fig. [5] and at 
that scale they do seem to coincide with points where Tq 
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FIG. 14: (Color online) Tail exponents Xq for eigenvectors 
of the random intermediate map for 7 — 1/3 with A'' = 2^^ 
(blue filled circles) and A'^ = 2^* (blue empty circles) and for 
7 = 1/11 with A'^ — 2^* (red squares). Same parameters as in 
Fig. O Inset: Tail exponent for 7 = 1/11 for larger values of 
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FIG. 15: (Color online) Tail exponents Xq for iterates of the 
Anderson map. Same parameters as in Fig. U The inset 
shows that 1/xq is well-fitted by a parabola. 



Xq = 1 is reached around g w 3 while Fig. [10] seems to in- 
dicate a separation of the multifractal exponents around 
g w 1. We note that for the finite sizes considered, the 
value of Xq seems to become infinite as g — >■ 1 , indicating 
that in this regime the distribution of moments is not 
any more fitted by a power law at large moments (see 
also Fig. [TO)) . The behavior of the exponents Xq will be 
further discussed in the next Subsection. 



E. Relation between multifractal exponents and 
moment distribution 



As explained in Section [ITl it was proposed in [21| 
that the exponents Tq and t'^p were related to the mo- 
ment distribution through the relation ([6]). This for- 
mula was proved only in some very specific cases, such 
as the PRBM model in the regime of weak multifractal- 
ity, but conjectured to be generically valid. The results 
of the preceding subsections enable to check numerically 
whether this formula holds for our models. 



1.25 



Tqxq 0.75 



0.5 



0.25 



















^- 












^■^ 




«,'' 


- 






^■^ 




>'' 


■ 


,o' 








,<3 




^''O 




'' «• 






- 



0.5 



0.75 



1.25 



XqTq 



typ 



FIG. 16: (Color online) Relation between Tq^^ and XqTq^^ for 
eigenvectors of the intermediate map for 7 — 1/3. Parameter 
values are the same as in Fig. [3] Blue/green (black/grey) 
circles correspond to values of q larger/smaller than 2.5 (the 
values closer to g = 2.5 are on the left for both series of 
points). Dashed line is the formula ([6]). 



and T*yP separate. However these points are far beyond 
the value g « 1 at which the multifractal dimensions Dq 
and D^yP separate at the scale of Fig. [TOl and also differ- 
ent from the values obtained by fixing the relative differ- 
ence of the exponents in Fig. [TT] (equal to (7+0.01 = 1-31 
for 7 = 1/3 and q+o.oi = 2.51 for 7 = 1/11). The value 
where the curves separate is dominated by the rare events 
in the extreme tail of the distribution. The different val- 
ues obtained indicate that indeed the numerical results 
are still far from the asymptotic regime. Similar con- 
clusions can be drawn from Fig. 1151 which presents the 
power-law tail exponents Xq obtained for the moment 
distribution of the Anderson vectors: the point where 



In Fig. ll6|17l wc show Tg^^ as a function of XqT^^^ for 
the random intermediate map with parameters 7 = 1/3 
and 7 = 1/11. For 7 = 1/3 (Fig. [T6[) . it shows a certain 
agreement with the conjectured law for small values of 
q. Similarly for 7 = 1/11 (Fig. [TT]) the agreement with 
the law ^ is good, but the actual slope seems slightly 
different. On the other hand, for 7 = 1/3 and larger 
values of q the relation breaks down. We note that there 
is a certain ambiguity in the formula, since as can be seen 
in Fig. [TH] one can have two values of Tqx^ for the same 
value of XqTg'^^ (corresponding to two different values of 
q). The results indicate that the relation ^ can indeed 
be seen, even if approximately, in other systems that in 
Anderson transition models. Interestingly, the case 7 = 
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FIG. 17: (Color online) Relation between Tqx^ and XqT^^^ for 
eigenvectors of the intermediate map for 7 = 1/11. Parameter 
values are the same as in Fig. O The blue solid line shows a 
linear fit of the data {y — 1.143; — 0.56). Dashed line is the 
formula ©. 



FIG. 18: (Color online) Relation between Tq and r*^^ for iter- 
ates of the Anderson map with k = 1.81, iV = 2", \tp{0)) = \i) 
with i = 1024 and t = 10*. Parameter values are the same 
as in Fig. [l] The blue solid line correspond to y — x — 0.29. 
Dashed line is the formula ((6|. 



1/11 corresponds to a case of weak multifractality. This 
might indicate that the relation is better verified in the 
case of weak multifractality, and only approximate in the 
general case. But we cannot exclude that the regime of 
weak multifractality leads to weaker finite-size effects and 
that the results for 7 = 1/3 would eventually converge to 
the law (|6]) for larger sizes and many more realizations. 
An additional problem concerns the different scales of the 
two figures \W\ and 1171 As the multifractality is weaker 
in the case 7 = 1/11, the values of Xq and Tq are larger, 
leading to a much larger scale for the data in Fig. [iTl It 
is possible that the finite-size effects are comparable but 
show more markedly in Fig. [16] due to its much smaller 
scale. 

In Fig. [T5] we present the same numerical analysis for 
the Anderson map. The results show that a linear law 
similar to ([6]) can be seen. The slope is close to one, but 
the curve is shifted by a relatively large offset (w 0.3). 
Note that the formula © was predicted for eigenvectors 
of the Anderson model and PRBM; here we are look- 
ing at iterates of wave packets, which can show different 
behavior. 

Interesting properties of the exponents Xq can be de- 
duced from the relation ^. As mentioned in the in- 
troduction, if such a relation is verified it implies that 
for q > q^ (or equivalently a;^ < 1) the inverse of 
the exponents Xq should follow a linear law. Indeed, 
for q > q^ the exponent t'^p is linear and one has 
qxqa+ = Tqx^ = Dqj.^{qxq - 1), thus z = qxq is solu- 
tion of the equation Dz = j^o;+- If this equation has 
an unique solution, then the quantity qXq is a constant, 
equal to q^ for q ~ q^^-, and thus 1/xq should be linear 
as a function of q for q > q^. In order to check whether 
this holds in our case, we plot in Fig. ^\ the values of 
1/Xq as a function of q for the intermediate random map 
with parameters 7 = 1/3 and 7 = 1/11. In both cases. 



a linear law agrees well with the data at large q, but 
obeys an equation different from the one predicted, with 
in particular an extra constant term which depends on 

7- 

A second consequence of relation ([S]) arises in the case 
where Dq is given by a linear function. This is in par- 
ticular the case for the intermediate map in the regime 
of weak multifractality. In this regime it was observed 
[la | that for not too large q the multifractal exponents 
Dq are very closely given by the linear approximation 
Dq = I — q/b. Inserting this relation into ([6]), we get 
that for q G]g_,(7+[ (where Dq and D^^p are equal) the 
exponents are given by Xq ~ b/q^ provided qXq also be- 
longs to ]g_,g+[. We note that this in turn predicts a 
value g_|_ — vb for the separation point between the two 
sets of exponents, contrary to the linear scaling found in 
the data shown in Fig. [TlJ According to these considera- 
tions, in the small q regime and for weak multifractality, 
a quadratic behavior of 1/xq should be observed pro- 
vided the linear regime extends beyond the point q+. In 
the case 7 = l/H for the intermediate map, the linear 
regime is verified quite far away from zero but breaks 
down before q^ (see Fig. fTU]) . thus Dq^ is not a lin- 
ear function of its argument. Still, Fig. ^§\ (lower panel) 
shows that the parabolic behavior of l/xq is retained for 
small values of q; the inset shows that indeed a quadratic 
fit is much better than a linear fit in this range. In the 
strong multifractality regime, Dq is not linear any more 
beyond q ~ I; in that case, the inset of the top panel in 
Fig. [TOl shows that the behavior of l/xq is linear rather 
than quadratic. 
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FIG. 19: (Color online) Inverse of the tail exponents Xq for 
eigenvectors of the random intermediate map. Top panel: 
7 = 1/3; blue and black dashed lines are linear fits in two 
different q range; dashed line (linear fit for large q) is 1/xq — 
0.36g + 0.23; inset is a blow-up of the small q range. Bottom 
panel: 7 — 1/11; blue and black dashed lines are respectively 
quadratic and linear fits in two different q range; dashed line 
(linear fit for large q) is 1/Xq = 0.36g — 0.55; inset is a blow- 
up of the small q range, showing in red the best linear fit 
for the same data. In both cases, the two ends of the error 
bars correspond to two values of Xq obtained from a fit of the 
moments probability distribution over two different intervals. 
Parameter values are the same as in Fig. [31 



F. Symmetry between exponents 

As described in Section [III it was predicted analyti- 
cally and observed in the Anderson model at the tran- 
sition that a symmetry exists between multifractal ex- 
ponents. Indeed, the quantity A^ ^ Tg — q + 1 follows 
the law Ag = Ai_g (relation ([71)). This relation was 
predicted on very general grounds, and expected to hold 
for all multifractal quantum systems. It was seen in the 
PRBM model [H], the Anderson model and ultra- 
metric random matrices [lO[. It was also predicted to 
occur in simple multifractal cascade models in [23 • We 
performed systematic calculation of the quantity Ag for 
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FIG. 20: (Color online) Anomalous exponents A, (blue 
dashed curve) and Ai_g (red solid curve) for eigenvectors 
of the intermediate map with random phases, 7 = 1/3 and 
jY _ 2^*. Parameter values are the same as in Fig. [3l The 
black thin line shows the parabola g(l — g)/6 corresponding 
to the linear regime. 



Figure [201 shows the results of this analysis for the ran- 
dom intermediate map. The presence of a large linear 
regime complicates the picture, since the linear law de- 
scribed above in Subsection IV Bl verifies the symmetry. 
Thus the intermediate map can show deviations from 
the symmetry only outside the linear regime. It turns 
out (comparing Figs. [7] and [201 that the symmetry ^ 
is only present in the linear regime and does not extend 
beyond its validity. This seems to indicate that the sym- 
metry is absent from these models. 

As the linear regime is much smaller in the case of the 
intermediate map without random phase, if the symme- 
try does not hold in this system we should expect a larger 
discrepancy in the nonrandom case. Fig. [21] displays Ag 
and Ai„g for the nonrandom intermediate case, showing 
that indeed the symmetry is verified in an even smaller 
range of q values than for the random case, in agreement 
with the small extent of the linear regime. 

Finally, we have also computed the quantities Ag and 
Ai_g for iterates of the random intermediate map and 
of the random Anderson map. As said in Subsection 
V.A, the precision of the exponents degrades for g < — 2 
in the case of the Anderson map, so verification of the 
symmetry relation is delicate. Nevertheless, our results 
indicated that while for iterates of the random interme- 
diate map the symmetry ^ still does not hold beyond 
the linear regime, in the case of the Anderson map the 
symmetry remains valid within the numerical error bars 
(which are however quite large) well beyond the linear 
regime (data not shown). This seems to indicate further 
that the symmetry is a feature of the Anderson model, 
which is clearly absent from the intermediate map. 

Although we cannot prove rigorously that the symme- 
try relation ^\ does not hold for intermediate systems, 
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FIG. 21: (Color online) Anomalous exponents A, (blue 
dashed curve) and Ai_q (red solid curve) for eigenvectors 
of the intermediate map (|10p with 7 = 1/3 and A*' = 2^^*. 
The (light blue, light red) shaded region indicates standard 
error in the least squares fitting. The multifractal analysis 
was done with box sizes ranging from 16 to 1024. The black 
thin line shows the parabola q{\ — q)/b corresponding to the 
linear regime. 



our results strongly indicate that it is violated in these 
systems as soon as the linear regime breaks down. To 
further confirm that our numerical method is able to ob- 
serve the symmetry ([7]) in a system where it is present, we 
computed A^ and Ai_g for ultrametric random matrices 
where the relation is known to hold [10| . Our numerical 
method was able to confirm unambiguously the presence 
of the symmetry in this specific case (data not shown) . 

We note that in [g^ the presence of the symmetry ^ 
for the Anderson transition was theoretically predicted 
on the basis of a renormalization group flow whose limit 
corresponds to a nonlinear sigma model. It would be 
interesting to see if a different nonlinear sigma model 
can apply to the intermediate map, or if it is the sign 
that these models cannot describe certain aspects of these 
systems. 



VI. CONCLUSION 

In this paper, we have studied the different multifractal 
exponents one can extract from the wave functions of the 
intermediate map and the Anderson map. Both models 
are one-dimensional, and thus allow much larger system 
size than the 3D Anderson transition, but in contrast to 
Random Matrix models such as the PRBM they corre- 
spond to physical systems with an underlying dynamics. 

Our results enabled to extract the exponents over a 
large range of q values for the intermediate map. We have 
checked that two methods widely used in other contexts 



(classical multifractal systems) , namely the box-counting 
and the wavelet method, can be used to obtain reliably 
the exponents, giving similar results in most cases, al- 
though the box-counting method seems more robust for 
large negative values. 

Our numerical data allow to confirm that the Anderson 
map introduced inll9[ and experimentally implemented 
with cold atoms |17| indeed displays multifractal proper- 
ties at the transition point, although this multifractality 
is weak. As concerns the intermediate map, our data 
confirm the existence of a linear regime for the multi- 
fractal exponents Dq and D^J^^ which was first seen in 
[l6| , well beyond the regime of weak multifractality. In- 
terestingly enough, the linear regime is much smaller for 
the nonrandom intermediate map. We checked that the 
exponents Dq and D^J^ are different also in the case of 
the intermediate map, even in the nonrandom case where 
no disorder is present. 

Our numerical study of the moments of the wave 
functions and the multifractal exponents show that the 
generic behavior of Dq and £)*''p predicted for the Ander- 
son transition [6| are present for the Anderson map. Our 
results enable to extract the value of g+ and g_ through 
the behavior of the moments of the wave functions to- 
gether with the values of a+ and a_; the fact that the 
value is different from the one obtained by direct com- 
putation of the multifractal exponents shows that finite 
size effects persist in such systems up to very large sizes. 
Note, however, that as our numerical computations cor- 
respond to very large system sizes, this might indicate 
that the asymptotic limit may be difficult to reach even in 
experimental situations. In addition, our investigations 
show that the relation ([6]) between the moments and the 
exponents conjectured in |2l| is only approximately ver- 
ified in our systems, even in the Anderson map. At last, 
the exact symmetry relation ([7]) between the multifractal 
exponents of the Anderson transition discovered in [2^ 
is not present in intermediate systems. 

Our results indicate that intermediate systems, and 
more generally quantum pseudointegrable systems, rep- 
resent a type of model with some similarities with the 
Anderson transition model of condensed matter, but with 
specific properties. In particular, the absence of the sym- 
metry present in the Anderson model between the expo- 
nents suggests proceeding with care in using the nonlin- 
ear sigma model to predict properties of these systems. 

We think that further studies of these two kind of quan- 
tum systems are needed in order to elucidate the multi- 
fractal properties of quantum systems and their link with 
spectral statistics. 
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